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A quantum mechanical superposition of a long-lived, localized plionon and a matter excitation 
is described. We identify a realization in strained silicon: a low-lying donor transition (P or Li) 
driven solely by acoustic phonons at wavelengths where high-Q phonon cavities can be built. This 
phonon-matter resonance is shown to enter the strongly-coupled regime where the "vacuum" Rabi 
frequency exceeds the spontaneous phonon emission into non-cavity modes, phonon leakage from 
the cavity, and phonon anharmonicity and scattering. We introduce a micropillar distributed Bragg 
reflector Si/Ge cavity, where Q ~ 10 5 — 10 6 and mode volumes V < 25A 3 are reachable. These results 
indicate that single or many-body devices based on these systems are experimentally realizable. 



Cavity-quantum electrodynamics (cQED) refers to the 
interaction of a single-mode of the electromagnetic field 
with a dipole emitter. cQED has provided new ways 
of controlling photons and matter (atoms, qubits, etc.) 
in both atomic and solid-state systems. The pro- 
gression from an atom "dressed" with a cavity pho- 
ton ("traditional" cavity-QED) [TJ |5] to a semiconductor 
microcavity-polariton (exciton plus cavity-photon) [3] to 
solid-state many-body polaritonic devices [H El E] has 
opened up new avenues for physical investigation as well 
as technology (e.g., single photon sources, novel lasers, 
long-range entanglement, quantum simulation). Moti- 
vated by this, we seek an analogous progression utilizing 
quantum sound instead of light. 

Phonons are more suitable for some tasks than pho- 
tons due to their slower speeds and smaller wavelengths 
(e.g., in signal processing, sensing, or nanoscale imag- 
ing). Our work builds off recent experimental results in 
nano-optomechanical systems [7], where cooling, coher- 
ent control, and lasing of mechanical vibrations has been 
achieved, as well as previous consideration of phonons 
as decoherence pathways [5], as tools for coupling quan- 
tum systems [SHU]) and even as a means for simulating 
many-body dynamics |12| . 

In this Letter, we show that a phonon-based analogue 
of the cavity-polariton is possible. Introducing a suitable 
high-Q phonon cavity, we calculate the cavity phonon 
coupling to a two-level system (TLS) in silicon (Fig. 1) 
and also losses due to spontaneous phonon emission from 
the donor into non-cavity modes, phonon leakage from 
the cavity, and phonon anharmonicity and scattering. 
Despite the phonon's dependent nature on its host mate- 
rial and the different (non-dipole) donor/phonon interac- 
tion, a strong coupling regime can be established, similar 
to cQED, where the phonon-TLS states are hybridized. 
The result of this mixing of cavity-phonon and matter 
excitation we term the c&vity-phoniton |13j . 

Implementation. Silicon is a promising candidate for 
constructing a cavity-phoniton system. The physics of 
shallow donors in Si have been understood since 1950s 
and experimentally verified, while transitions between 
low energy donor states are known to be driven by acous- 
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Figure 1. (a) A cavity-phoniton can be constructed in 
a Si/Ge heterostructure cavity as a hybridized state of a 
trapped single phonon mode and a donor TLS placed at a 
maxima of the phonon field, (b) The P:Si donor lowest Is 
valley states, Ai, T2, and upper levels (Is/ 2p); their energy 
splittings can be controlled by the applied strain in the Si 
cavity, (c) Angular dependence of the coupling g q (0), Eq.|2j, 
for the deformation potentials of Ref. |15| vs. dipole ~ cos#- 
dependence (thin circles). 

tic phonons [5]. The six-fold degeneracy due to Si's 
multi-valley conduction band is lifted both by applied 
strain (e.g., due to the lattice mismatch with a substrate) 
and the sharp donor potential. Crucially, in [001] com- 
pressively strained-Si (Fig. 1), the first excited state at 
zero magnetic field of a phosphorous donor approaches 
~ 3.02 meV (0.73 THz): a so-called excited "val- 
ley" state. (The excited valley state has an s-like en- 
velope function like the ground state but opposite par- 
ity; because of this, valley state relaxation times can 
be much longer than for charge states). The energy 
splitting implies longitudinal (transverse) wavelengths of 
A; ~ 12.3 nm (At ~ 7.4 nm). For comparison, the en- 
ergy splitting to the upper 2p-like state is > 30meV 
(A2 P ~ 1.2 nm), unlikely to be amenable to phonon cav- 
ities. Since the P:Si Bohr radius is a* B < 2.5 nm in the 
bulk, A > a* B allows for easier donor placement, avoid- 
ance of interface physics, and bulk-like wave functions. 

A prototype implementation of a cavity-phoniton sys- 
tem in silicon is sketched in Fig. la. A strained-Si phonon 
cavity grown in the [001] direction, of length d c ~ A and 
lateral size D, is enclosed by acoustic DBRs formed as 
layered, epitaxially-grown, and strain-relaxed SiGe super 
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lattice (SL) heterostructures. The cavity length is cho- 
sen to be less than the critical thickness due to strain 
(see, e.g., Ref. [10]). The DBR SL unit periods con- 
sist of subsequent layers (A, B) — Si\- XAB Ge XAB of 
thickness (d^de), strain matched to a Sii_ s Ge s sub- 
strate (xa( B j = 0.55(0.05), s = 0.26, maximizing con- 
finement) and an appropriate capping layer, depending 
on the actually confined phonon mode. Note that ID 
DBR SL phonon cavities (D ^> d c ) are well understood 
and have been demonstrated in THz phonon cavities in 
III-Vs O [T3]; coherent phonons in SiGe superlattices 
were studied as well [20J. In the case of micropillar DBR 
(mpDBR) structures (designed to increase the phonon- 
donor coupling), the DBR lateral dimension may become 
comparable to the phonon wavelength {D > d c ~ A) 
and the confined mode is a mixed longitudinal/ transverse 
one. The trapped mode with wavelength A g and phase 
velocity v q is designed to be resonant in energy with the 
first excited state of the donor. Similar to the ID DBR 
[5] , the thickness of the SL unit cell is set to match the 
Bragg condition d 
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4 B A„/4v„, where v$ B are the 
phase velocities (using isotropic approximation, see, e.g., 
Ref. [13]). The donor is placed at the center of a A-cavity 
{d c = Ag), where the displacement u(r) is maximal [24 . 

Hamiltonian and coupling. In the semi-classical pic- 
ture an acoustic phonon creates a time-dependent strain, 

£ a p(r) = \ (f^ + which modulates the energy 

bands and can drive transitions in a localized state, e.g., 
a donor. For Si, from the multi valley electron- phonon 
interaction US] one can derive the matrix element be- 
tween valley states, \s,j): 
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where k\j are the directions toward the valleys, s, s' 
label the orbital (envelope) function(s), and H„(Si) ~ 
8.77 eV, Srf(Si) « 5eV [T5] are deformation potential 
constants. For the donor-phonon Hamiltonian we ob- 
tain the interaction (of Jaynes-Cummings type): H g ss 
hg q {o-+ e b qa + cr~ e b q a ) where only the resonant cavity 
phonon with quantum numbers q, a and energy hw q a is 
retained, b q a is the phonon creation operator, and <7+ = 
\e){g\ refers to the donor transition between ground and 
excited states. In the loss part: H\ oss = H K + H' a , + Hr, 
H K couples the cavity mode to external continuum of 
other modes giving a cavity decay rate k — Lo q:<y /Q (ex- 
pressed through the Q- factor); H' anh includes phonon de- 
cay due to phonon self-interaction and also phonon scat- 
tering off impurities (mainly mass fluctuations in natural 
Si) . The coupling of the donor to modes other than the 
cavity mode, Hr, leads to its spontaneous decay. 

The valley states: ls(A 1 ), ls(T 2 ), that make up the 
TLS are the symmetric and anti-symmetric combina- 
tions of the conduction band valley minima along the 



z-direction (see Fig. lb). Due to opposite parity of the 
states the intravalley contributions cancel. The inter- 
valley transitions are preferentially driven by Umklapp 
phonons [3] with a wave vector q at q u ~ 0.3 where 
q u = G+i — 2k z is the wave vector "deficiency" of the 
intervalley k z — > —k^ transition, and G+i = ^(0,0,1) 
is the reciprocal vector along z. Since typical values give 
q u r s» q u a* B ~ 9.4 > 1 and qr ~ 1 for 3meV, the coupling 
is calculated exactly (not using the dipole approximation) 
for longitudinal and transverse polarizations 
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where a G w 0.3, I ge {9) = J dr[<J>f s (r)] 2 e -"f sin (q u r) = 
— 2 fiq cosQ (i-7 g cos g) — . intervalley overlapping 

factor, a q = 1 + \{q 2 a 2 + cfcbr 2 ), f3 q = ^b 2 qq u , ~/ q = 

j^-(a 2 — b 2 )q 2 , and a/b are the radii of the Kohn- 

Luttinger envelope function <&f s (r) (see, e.g. Ref.|15|L 
The calculated coupling is to plane wave modes, related 
to a rectangular cavity with periodic boundary condi- 
tions. Fig. lc shows the directionality of the coupling for 
longitudinal and transverse phonons. The angular de- 
pendence in Eq. ^ for longitudinal phonons is similar 
to dipole emission (I^ p ~ cos 9) , but enhanced in a cone 
around indirection due to non-dipole contributions. Un- 
certainties in S<j calculations (see Ref. [5S]) result in an 
overall factor of 2 difference in the maximal coupling. 

Coupling to cavity mode. The matrix element hgq is 
just the interaction energy of the donor with the "phonon 
vacuum field" and expresses a generic dependence oc 
1/VV on the normalization volume; it goes to zero for 
large volumes. In a cavity, V is the physical volume of 
the mode [3]. By virtue of Eq. we first consider a 
DBR cavity with a length d c — X\ ~ 12.3 nm designed 
for longitudinal resonant phonon along the z-direction 
(we use isotropic velocities for Si: vi — 8.99 10 3 ms _1 , 
Vt = 5.4 10 3 m s -1 , see Ref.[23j). Taking the minimal 
lateral size D m i n ~ A/ to ensure D > 2a* B rj 5 nm, one 
gets a mode volume of V m i n s» A 3 for P:Si. Thus, we 
estimate the maximal phonon-donor coupling as gw — 
3.7 10 9 s _1 . For D — 5 A the coupling is still appre- 
ciable: = 7.410 8 s _1 . Surface undulation typical 
of step-graded SiGe quantum wells (see Ref. [T^) gives 
D w 200 nm < 15A, though this interface imperfection 
is avoidable with heterostructures grown on defect-free 
nanomcmbranc substrates [26J. 

For a realistic cylindrical mpDBR cavity, the modes 
can be constructed as standing waves with energy Hui, 
and wave number q along the pillar z-direction, with 
stress- free boundary conditions on the cylindrical surface. 
The displacements for compressional modes (see e.g. 
Ref. jllp are: u r [r,z) = [A r Jifar) + B r Ji(rjtr)] sinqz, 
u z (r,z) = [A z J a (r]ir) + B z J (n t r)} cos qz, whe re J ,i(r) 

are Bessels of 1st kind, rjij = ^Juj 2 /vf t — q 2 , and A t 
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Bi are constants. These modes are lower in energy and 
couple strongly to the donor; the related strain has a 
node at the Si-cavity z-boundaries (for A-cavity) |24| . 
For a fixed resonant frequency u) and lateral size D, the 
dispersion relation q — uj/v q (D) has multiple solutions 
qi,i = 0,1,2,..., where qo stands for the fundamental 
mode, q\ for the 1st excited mode, etc. Each mode prop- 
agates with its own phase velocity v qi {D) ^ vi,vt- For 
a A-cavity and D = A/ we calculate via Eq.Q maximal 
coupling to the fundamental mode with X qa — 6.9 nm 
to be g$ — (7 max = 6.5 10 9 s _1 (see Table I), com- 
parable to the above estimation. For larger D 7 how- 
ever, the coupling to the fundamental mode rapidly de- 
creases (e.g, gf^ = 5.8 10 5 s~ 1 ) since the mode trans- 
forms to a surface-like Rayleigh wave (vR ay ieigh < v t ). 
Coupling to higher mode branches is appreciable and de- 
creases roughly as l/\/V, e.g., for the 1st mode branch, 
g$ = 2.4 10 9 s" 1 and g$ = 3 10 8 s~ 1 . Among vari- 
ous mode choices we note that for any diameter D there 
is a higher excited mode with resonant wavelength close 
to that for longitudinal phonons. E.g., for D = 3A/ we 
found the wavelength of the resonant 4 th excited mode 
as X q4 ~ 12.4nm and the coupling g^ = 3 10 9 s _1 (i.e., 
values similar to the rectangular DBR estimate). 

Loss. Losses in this system are dominated by donor re- 
laxation and leakage of the confined phonon mode. Sim- 
ilar to cQED (see, e.g., Ref. [3]), one can argue that 
the donor relaxation, r rc i ax , to modes different than the 
cavity mode (and generally not trapped into the cavity) 
is bounded by the donor spontaneous emission rate in 
the bulk: r re i ax < V. We calculated the bulk donor re- 
laxation to longitudinal and transverse phonons to be 
(0 - nn7 ,-i a „A r(*) - 9 2 10 7 s' 1 , respectively, 

(0 i rW _ioin8 s -ly 



Y ( gl = 3 10 7 s- 1 and T ge 



T g l = 1.2 10 s 



for the 3meV transition (r = T ge 
The relaxation to photons is electric dipole forbidden and 
suppressed (Hj by (\ pho ton / a* B ) 2 ~ 10 10 - 

The cavity mode loss rate is calculated as mainly due 
to leakage through the DBR mirrors, similar to optical 
DBR cavites [55] (except that for phonons there is no 
leakage through the sides). Generally, for the cylindri- 
cal micropillar DBRs, the cavity mode involves coupled 
propagation along the micropillar of two displacement 
components, u z (r,z), u r (r,z), and two stress fields, T zz , 
T zr . For small diameters, D <C A;, the fundamental mode 
becomes mainly longitudinal (~ u z ), propagating with 
the Young velocity, vq = \JEjp. Using 4x4 transfer 
matrices we calculated QmpDBR — 10 6 for N = 33 layers 
for the confined, mixed fundamental mode at D = Af, 
that is close to the limiting Q-factor related to pure lon- 
gitudinal propagation|31 . (This is also similar to the ID 
DBR value relevant for D 3> A;). For our design we ob- 
tain a cavity loss rate k = A v /hQ ~ 2.8 10 6 s _1 . This 
can be decreased by adding more layers. 

At low temperatures the phonon anharmonicity losses 
are negligible (a rate r an h — 1.4 10 4 s _1 at 3meV), 
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r 2 /2 3 2 


< 10" 6 


< 210" 4 


< 610~ 7 


# Rabi flops 


2 S /(«+r) 


~ 100 


~ 102 


~ 93 



Table I. Key rates and parameters for circuit QED |33j (ID 
cavity) vs. the phoniton system; we show calculations for 
maximal coupling for a A-cavity with lateral diameter of D = 
A; (in general, g ~ Q — 10 6 , and comparable number 

of Rabi flops |37| . 



while scattering off impurity mass fluctuations in nat- 
ural Si amounts to a rate two orders of magnitude larger: 



imp 



710 s . It is notable that in isotopically pu- 
rified bulk silicon (an enrichment of 28 Si to 99%) the 
scattering rate will decrease by an order of magnitude 
and the related phonon mean free path will be of the 
order of vi/T- lmp ps 10 cm |11| l29l [32] . In this case, 
the cavity leakage dominates, k 3> {r an h,r imp }, and 
the number of vacuum Rabi flops can reach as high as 
"-Rabi = 2fj(D)/(r + k) ~ 102 for a cavity Q-factor, 
Q = 10 6 , and some n Rabi ~ 77 for Q = 10 5 . For 
D = 10A; and similar Q one still has: n\ll hi — 1 (17) for 
the 1st (2nd) excited mode. Further, nearby modes can 
be well separated from the resonant mode, e.g. for the 
fundamental mode and D = A; the next mode (in trans- 
verse direction) is ~ 0.3 = 0.9 meV off; the transverse 
separation for D = 10A; (for the 1st excited mode) gives 
0.009 A v , which is more than two orders of magnitude 
larger than the linewidth To ~ (r + k)/2 of the two hy- 
bridized levels. 

Experiment and Discussion. We have shown that the 
donor:Si cavity-phoniton can enter the strong coupling 
regime with 2g/(T + k) ~ 10 - 100 (Table [[]). A principle 
experimental confirmation would be the observation of 

the "vacuum" Rabi splitting: fl = [g 2 - (T - k) 2 /4] 1,2 ; 
two resolved spectral peaks can be observed if 2Q > T . 
The Rabi splitting can be enhanced as f2 ~ gVN by 
placing more than one donor (N > 1) in the cavity [TJ |3] 
(e.g., via a delta-doped layer). This could allow for large 
coupling even for large diameter micropillars/ ID DBR 
structures (since ^£ cx = 1). Further, strain or 

electric field (from a top gate) [TTJl IMl 155] can be used to 
tune the valley transition into resonance. 

Experimental techniques are available to probe the Si- 
phoniton (low temperature, T ~ IK, and low phonon 
numbers are assumed). First, free-electron lasers have 
been used to probe the Is — 2p transitions in P:Si |18| . 
Observation of the vacuum Rabi splitting is possible by 
measurement of the absorption spectrum of the allowed 
optically probed transition ls(l2) — > 2po (~ 30meV) us- 
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ing weak optical excitation. Appropriate phonons can be 
introduced to the system by excited valley state emission, 
by piezo-actuators, or by increasing the temperature, as 
was done for the first observations of Rabi oscillations [36 
(phonons of 3 meV ~ 30 K). Second, pump and probe 
optical techniques have been demonstrated to observe co- 
herent phonon effects in III-V [7\ US] and SiGe |2D] SL 
heterostructures. Observing the reflected phonons from 
this structure will show the phonon-Rabi splitting char- 
acteristic of cavity-QED systems pQ . 

The cavity-phoniton can be realized in other materials 
and systems. In particular, our system should be compat- 
ible with recently demonstrated (though in the few GHz 
range) high-Q phononic band-gap nano-mechanical and 
opto-mechanical (membrane) cavities in silicon (e.g., in 
[7] and [35]) ■ Quantum dots, spin transitions, color cen- 
ters in diamond, and other donors (particularly Li:Si [10 ) 
may offer smaller resonance energies and correspondingly 
larger cavities (wavelengths) [32] ■ In the case of [001]- 
strained Si considered in this paper, the two lowest lev- 
els in Li possess essentially the same state structure as 
P:Si and approach a splitting of Au — 0.586 meV for 
high strain (from zero splitting at no strain). We calcu- 
late (see Table |l| the corresponding Li:Si donor-phonon 
coupling for the D = A; (now A; = 63.2 nm) refer- 
ence cavity; strong coupling can still be reached. For 
the DBR cavities, the Si critical thickness can be made 
80 — 100 nm by lowering the Ge content in the substrate. 
For 2D-phononic bandgap cavities, direct numerical cal- 
culations for cavity-trapped phonons (u> r /2ir ss 10 GHz) 
in novel Si-nanostructures [35] show the potential to 
reach Q cav > 10 7 in the ideal case. 

The phoniton is a new component for constructing 
and controlling macroscopic artificial quantum systems 
based on sound. Besides single phonon devices, sys- 
tems composed of many coupled phonitons could ex- 
hibit novel quantum many-body behavior. For exam- 
ple, "solid-sound" systems in analogy with coupled cavity- 
QED solid-light systems [5 could demonstrate Mott in- 
sulator like states of phonons in coupled phoniton cavi- 
ties. Cavity/qubit geometries such as these may also be 
relevant for quantum computing (QC) applications: to 
mediate interactions between distant qubits or inhibit de- 
coherence. The systems proposed here will benefit from 
the drive in silicon QC towards more purified materials, 
perfect surfaces, and precise donor placement. 
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SUPPLEMENTAL MATERIAL 
Deriving the phoniton Hamiltonian 

In semiconductors at low energies, the conduction band 
plays a role similar to the QED vacuum: the quasi-free 
electrons propagate according to a quadratic dispersion, 
scatter off acoustic phonons that possess linear disper- 
sion, and can be bound to charged impurity potentials 
to establish hydrogen- like atoms (donors). Thus, the 
pattern of QED (matter and photon interactions) re- 
peats itself in the behavior of electronic quasiparticles 
in the solid (the "matter") and the phonon vibrations in 
a crystal (the "photons"). Consider the electron-phonon 
interaction. The Hamiltonian can be rearranged into the 
electron part H c , the phonon part -ff p h, and the electron- 
phonon interaction i? .ph: H = H c + _ff p h + -ffc.ph, where 

2 

JJ C = + V\ at {r) includes the periodic lattice poten- 
tial Viatic) °f the perfect crystal acting on an electron at 
point r. For small atomic displacements it(r) an expan- 
sion in normal modes gives 



u(r) 



E 

q.a 



u qa (r)b qa + u* qa (r)b qa ) 



(SI) 



which approximately diagonalizes the phonon part H^: 
H ph = E, !0 H<t ( b q <j b q<? + \) + #anh- The small an- 
harmonicity, H an h — cb qa b q > a bk a +•••, is related to 
phonon self-interaction. The mode normalization in Eq. 
SI is f d 3 r u*„(r)u„„(r) = 7^ — , so that b$,„ creates 
a phonon in the mode q, a with energy huj qa in a ma- 
terial with mass density p. The vector q denotes a col- 
lective index of the discrete phonon mode defined via 
the phonon cavity boundary conditions and quantization 
volume V. In particular, the plane wave expansion corre- 
sponds to rectangular periodic boundary conditions and 
1/2 

zq r with wave vector q, po- 



V( r ) =" \2pV^„ ) *W 

larization £q,o-, and phonon branch a. 

By considering low-energy acoustic phonons for elec- 
trons close to the band minimum, k fa kg, the electron- 
phonon interaction can be written as 



(S2) 



The operator Djj — —pipj/m+Vij(r) (withp = —(i/H)V 
and a crystal model dependent Vy(r)) coincides with 
the constant-strain deformation potential [TJ [5] and the 

causes transitions between 



strain e, 



.( r \ = l(dui, duA 
IV' ) — 2 \dr-i + dr, J 



electronic states. A donor impurity center makes a per- 
turbation to the perfect crystal, thus binding the low- 
energy electrons. One introduces electron field oper- 
ators vp(r) = J2 S c s ips( r ) related to the donor bound 
states ip s (r), where c s is the annihilation operator for 



density operator, vf't (r)^(r), the electron-phonon in- 
teraction takes the second-quantized form: -ff C! ph = 

£*,,' I2q,a c l' c s( b q,<? + b q,<?) V qJ where the matrix el- 
ement, V^'j = (ip s >;{q,<j}\H™ ph \ip s ), describes phonon 
transitions between the bound states. 

Similar to quantum optics, the electron-phonon inter- 
action -ffe.ph describes transitions in an n-level atom (al- 
though the analog of electric-dipole approximation is not 
made). The Jaynes-Cummings type Hamiltonian [5] re- 
sults when a single phonon mode q is in resonance with 



the transition (i.e. 



J q,(T 



/). This is possible if we 



have a suitable phononic cavity and the phononic band- 
width Auj q is much smaller than the transition frequency 



Introducing rising (lowering) operators 



c\,c s 



(cr = c\c s i) for the specified transition (s — > s') and 
supposing real matrix elements such that V q § = V q s a = 
hg q , we obtain H g = hg q (a + + a~) (b qtlJ +b q cr ). In 
the typical case of g q <C oj ss i only the "energy con- 
serving" operators a + b qyCr , o~b qa will survive in the 
rotating wave approximation obtaining H a ss iJ JC — 



hg q 



a b qa ). The total Hamiltonian in such 



iq \y ^q,a 1 " "q !( 

a cavity approximation can then be written as 



H = hw. 



1 

2 



HA 



(S3) 



that state. By integrating Eq. ( S2 1 over the charge 



where we have denoted the atom transition frequency as 
A = uj ss < . 



Phonon-donor coupling 

For shallow donors the wave function is localized in 
"valleys" in fc-space, near the local conduction band min- 
ima kj, and for Si (tetrahedral symmetry) these are dis- 
placed at k ~ 0.85 27r/a along the equivalent directions 
x, —x, y, —y, z, —z in the qubic lattice. The valley 
state is a modulated Bloch wave: ip{(r) ss ipk\r)& s (r) 
with the well known envelope functions & s {r) Each 
donor bound state is superpositions of such valley states, 
ijj s { r ) = £j=i a j^ii r ) with a s j being the "valley pop- 
ulations". The six-fold degeneracy of conduction elec- 
trons in silicon is broken both by strain (e.g., due to 
the lattice mismatch with a Sii-^Ge^ substrate) and 
the sharp potential of the phosphorous (or other) donor. 
Thus the two lowest states in our device approach sym- 
metric and antisymmetric combinations of the ±z val- 
ley states (see Fig. 1, main text) [S] and asymptote to 
E s > s = E e - E g ~ 3.02 meV splitting for large (x > 0.1) 
mismatch. Both states are "s-like", having the same s- 
like envelope function but opposite parity due to differ- 
ent content of the highly oscillating Bloch components. 
The next highest state is "p-like" and lies some 30 meV 
above the ground state. 



7 



One can calculate the partial matrix element for inter- 
valley transitions in the long wave-length (acoustic) limit 
OH: 

1/2 



xJ2 a G I dr [®i>( r )Y & s (r) e~ l ^ +q + G - k ^ r . (S4) 

G ■* 

The total matrix element is given by a sum over all valley 
contributions weighted by the valley populations: = 
a t <Xj Vfj s (q, a). Here, E d and E u are the low-q de- 
formation potentials (that compose in Eq. S2 1 , fej is a 
unit vector towards the ith valley (e.g. k\ = x, &2 = —x, 
etc.), and are the Fourier series expansion coefficients 
of the function u* k , {r)uk j (r) (uk are the periodic parts 

of the Bloch waves), i.e. u^. (r) Ukj (r) = J2g a G e ~ lGr ■ 
In the sum over G the main contribution comes from 
the minimal reciprocal vectors. Higher terms are sup- 
pressed due to spatial confinement of the orbital part 
and/or rapid decreasing of the Fourier coefficients for 
higher G [7]. Due to different parity of the valley pop- 
ulations under inversion kg — > the intravalley con- 
tributions (i = j) will cancel. Leading contribution to 
Vq§ comes from intervalley (Umklapp) transitions with 
q x, q u = G+i - 2kg, where G ±1 = ^(0,0, ±1) is the 
reciprocal vector along the z-direction [&]. 

Donor relaxation 



The total emission rate is calculated from Eq. (S4) 



via the Golden Rule: T^J = £ q ^\V^\ 2 8{u q ^ - oj ge ). 

Using the notations x = cos#, J(x;f3 q ,j q ) = j^-I s °(x) 

(cf. Eq. 2 of the main text), and S^'(x) = + 5 u x 2 , 
S^(x) = E u xvi — x 2 we obtain the emission rates for 
acoustic modes: 



9C 47T phv 7 a 



1 + \ (? 2 « 2 + 
dx [~S°\x)J(x-f5 qnq )] 2 



(S5) 



where v a is the speed and qfe — U3 ge /v a is the wave 
number of the emitted phonon of polarization a. We 
mention that T q a J as a function of energy of the tran- 
sition must experience a maximum when qfe is close 



0.3 



This corre- 



to the Umklapp value q u or 
sponds to transition energies 19.2 meV and 13.3 meV for 
longitudinal and transverse phonons respectively, which 



is in accordance with Eq. ( S5 ) . Using the mass density 



p = 2330 kg m~ 3 for Si and the deformational poten- 
tial constants: S u (Si) ~ 8.77 eV, S^(Si) ~ 5.0 eV, and 
ac ~ 0.3 [7] one obtains the rates at 3meV: 

r« = 3 10 7 a" 1 , = 9.2 10 7 s" 1 . (S6) 

For reference, the values calculated in the dipole approx- 
imation are about w 35 — 50% larger than the above. 



Cavity design and cavity quality factors 

ID DBR cavity 

Due to the analogy between light and sound in solids, 
the theory of optical DBRs can be applied to the 
design of acoustic DBRs. For acoustic waves prop- 
agating perpendicularly to the ID DBR SL mirror, 
Bragg type interference leads to forming of frequency 
stop-bands |8J in the reflectivity spectrum, positioned 
around a central frequency (for the m th -order band) 
u) m = mir (cIa/va + cIb/vb) ■ The presence of a 
resonant Si-cavity (a "defect") in the SL leads to a 
dip in the reflectivity spectrum (cf. Fig. |2|, whose 
width would represent the ID DBR Q-factor, provid- 
ing the cavity losses are mainly through the DBR mir- 
rors. Similar to optics, the transfer matrix method (see, 
e.g., Ref. [9]) allows us to evaluate Qdbr — 2it{d c + 
i_DSfl)V^/Alni? c , where L DBR = Z A Z B \ q /2Z c \Z B - 
Za\ is the acoustic effective DBR mirror length and 
R c = [(Z S /Z C -Z 2N )/(Z S /Z C + Z 2N )] 2 is the DBR 
peak power reflectivity calculated for N layers with 
layer impedance ratio Z r — pava/ Pbvb {Z s /Z c are the 
impedances of the substrate/cavity respectively). 

The cavity quality factor extracted numerically from 
the linewidth of the reflectivity dip (see Fig. |2| agrees 
well with the simple analytics. For our design, we choose 
the layer materials A = Sio.45Geo.55 and B = Sio.95Geo.05 
with thicknesses of d A — 2.1 nm and d B — 2.8 nm (strain 
matched [10] to a substrate of Sio.74Geo.26), to confine 
a longitudinal phonon (d c = A;); we obtain a Q-factor 
of - 1.2 10 5 for N = 29 SL layers; for TV = 36 Qdbr 
reaches 10 6 . 



Micro-pillar DBR cavity 

While longitudinal and transverse phonons formally 
decouple under isotropic approximation in the bulk, 
they combine in mpDBR cavity/waveguide modes due 
to the strong effect of the side boundaries, and prop- 
agate with a differing phase velocity v q ^ vi,v t . For 
the cylindrical mpDBR cavity we focus on compressional 
modes only, since flexural modes are energetically higher, 
while torsional (transverse) modes do not couple to the 
donor. The boundary condition on the mpDBR cavity 
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36 layers : Q a 10 6 
Sio.45Geo.55 // Sio.95Geo.05 
d c = X = 11.6 nm 






0.7 

f (THz) 



-1. 0. 1. 




Figure 2. Reflectivity and field distribution of acoustic distributed Bragg reflector (DBR) phonon cavity, (left) Reflectivity 
spectrum of the phonon cavity enclosed by DBRs consisting of TV = 36 subsequent layers of thicknesses cLa — 2.1nm/ds ~ 2.8nrn 
(A = Sio.45Geo.55/-B = Sio.95Geo.05)- Slight thickness changes due to finite interlayer spacing approximation do not cause 
significant change in Q. For the ID calculations we use the [001] phonon velocity of vi — 8.4310 3 m/s giving A; ~ 11.6 nm. 
The sharp dip in the reflectivity at fo = 0.73 THz corresponds to the donor transition energy ~ 3meV. The left inset shows a 
zoom in of the reflectivity dip and corresponds to the confined cavity mode with FWHM of 0.6 MHz and Q w 10 6 . (middle) 
Amplitude of the localized phonon displacement in the ID DBR cavity is shown with respect to the position along [001] in 
Si, with the Si cavity of thickness d c and penetration depth S p into the DBR, respectively, (right) The mixed phonon mode 
in a 3D micro-pillar DBR structure. Normalized amplitudes of the phonon displacements u z (r, z) and u r (r, z) are shown for 
the fundamental mode inside the Si cavity of thickness 6.9 nm and diameter D — 12.3 nm. The DBR mirrors enclosing the 
Si cavity are constructed from alternating layers of Sio.95Geo.05/Sio.45Geo.55 with thicknesses of 1.7/1.2 nm, respectively. 3D 
cavity reflectivity spectrums are qualitatively similar to the ID case. 



z-interfaces is linked to the Bragg resonance condition at 
the first zone center stop band (for A g -cavity), leading 
to q = 2ir/d c = A v /Hv q , where q is the cavity phonon 
wave vector along the pillar z-axis. Additionally, the 
free-standing (stress-free) boundary condition along the 
cylindrical surface implies for the stresses: T rr = T rz = 
at r = R, and leads to the well-known Pochhammer- 
Chree dispersion relation of compressional waves (see, 
e.g., Ref. mi), 2r,i{rf t + q^Mm^MvtR)/ R ~ (v? ~ 
q 2 ) 2 J (mR)Ji(VtR) - 4q 2 ViVtJi(mR)Jo(VtR) = 0, where 
Jq.i(t) are Bessels of 1st kind, R = D/2 is the radius, 
and T)i = ^u) 2 /vf — q 2 , X] t — y/u 2 Jvf — q 2 may be inter- 
preted as transverse wave vector components. For fixed 
u, D, the dispersion relation possesses multiple solutions 
qi — uj/v qi (D), i — 1,2,.... Each mode branch defines 
its own phase velocity v qi , v qo < v qi < . . . with the fun- 
damental mode branch having the smallest velocity. 

The modes in the mpDBR cavity can be constructed 
as standing waves using waveguide solutions for compres- 
sional modes: 

2q 2 J^R) 



u r (r) = -Arji 



Ji{m r ) 



Mvtr) 



u z (r) = Aq 



Joivir) 



V 2 - q 2 M Vt R) 

Zmvt JiiviR) T , 

-2 ~2 j 7Z. m MVtr 



sm qz 



cos qz. 



The relevant strain components have ~ sin qz depen- 
dence, i.e., the strain has a node at the Si-cavity z- 
boundaries to ensure non-zero coupling with the donor 
A\ <H- T2 transition (see Eq. 2 of the main text) . 



The constant A can be calculated from a proper nor- 
malization condition: 



V ca , 



d 3 r u* q (r)u q (r) 
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2puj q 



(S7) 



where u q (r) = e r u r (r, z) + e z u z (r, z) is the displacement 



"vacuum field" (see Eq. SI I in the Si-cavity. Here 7 is 



the fraction of the phonon vacuum energy, Huj q /2, stored 
in the cavity of length d Cl while an energy fraction 1 — 7 
is stored in the DBR mirrors, related to the penetration 
dept S p shown on Fig. [2] Approximate transfer matrix 
calculations based on the propagation of four fields (the 
displacements u z , u r , and the stresses T zz , T zr ) along the 
mpDBR waveguide give 7 ps 0.3 (note that the coupling 
g scales as ~ ^/7 which renormalizes the g m ax's calcu- 
lated in the text) that is commensurate with the value 
estimated from the ID DBR field distribution. 

For each mode of mpDBR considered in resonance 
with the donor, a suitable cavity design with appropri- 
ate cavity design parameters, i.e. cavity and SL layer 
thicknesses, strain matched substrate Ge content, etc. 
is required. In accordance with our ID DBR design, 
layer materials A and B are kept the same. For a di- 
ameter of D = Xi = 12.3 nm, a Si-cavity thickness of 
d c = A 9o = 6.9 nm results with the fundamental mode be- 
ing in resonance with the P donor (the layer thicknesses 
are correspondingly = 1.7 nm and dg = 1.2 nm). Us- 
ing the TM method, the ideal Q-factor of the funda- 
mental mode is calculated to range from 10 4 — 10 6 for 
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TV = 20 — 33 unit layers of the DBR, and similar num- 
bers appear for the higher excited modes. 



Anharmonicity and scattering off impurities 

In isotropic crystals at low temperatures (T <C fiw/fce) 
a TA-phonon decay is forbidden, while the LA-phonon 
decays through LA -> TA + TA and LA ->■ TA + LA 
channels. For Si, r|°| « 4.5 10 4 i4 Hz M; this S ives for 3 
meV (0.73 THz) a rate 1.4 10 4 s" 1 (Si anisotropy lifts out 
the TA selection rule, leading to same order of magni- 
tude TA-decay rate; see, however, [TJ] ). For larger tem- 
peratures (T > 3meV/fcs ~ 30 K), Tla scales as wT 4 
(Landau-Rumer effect, see e.g., Ref. [11]) leading to an 
order of magnitude enhancement. Stronger phonon losses 
come from scattering from isotopic point defects. For 
natural Si the isotope mass fluctuations lead to a rate 

r imp w 2.4310 6 4hz C2. which S ives ~ 710 5 s~ 1 for 
0.73 THz. The estimated rate gives a phonon mean free 



path of the order of uz/Ti mp ~ 1.2 cm. An enrichment 
of 28 Si to 99% (natural abundance of 28 Si is 92.23%) de- 
creases the scattering rate by one order of magnitude; 
so, the dominating phonon loss may become scattering 
off the surfaces and interfaces of the phonon microcavity 

In Ref. [15 via direct measurement of the phonon 
decay rate in a GalnAs SL cavity {y\ — 0.46 THz): 
T p h — r con f + r an h, an anharmonicity rate was estimated 
as r an h — 3.310 9 s~ 4 , taking into account the phonon 
confinement due to the SL DBRs. The observed large 
rates are nevertheless in qualitative agreement with the 
above theoretical estimations, since the experiment was 



performed at a high (room) temperature and dealt with 
large number of initially injected phonons |11| . Thus, 
we expect that in quantum-limited experiments at low 
temperatures (~ 1 K) and low phonon occupation num- 
bers scattering rates of the order of 10 4 — 10 5 s~ 4 can be 
achieved. 
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